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Abstract 


A method for the optimization of structures to satisfy fiutter velocity 
constraints is presented along with a method for determining the flutter 
velocity. The material presented represents a summary of Ref. I through 5 
which are a direct result of the research effort of this grant. A method 
for the optimization of structures to satisfy divergence velocity constraints 
is presented in Appendix B. 



Acknowledgments 


The author wishes to thank Dr. Kumar G. Bhatia, Mr. Jerry L. Cooper, 
Mr. Salfee K. Motiwalla and Mr. Yee-Yeen Chu for their diligent assistance. 

The author also is indebted to Mr. Robert V. Doggett, Jr. of NASA 
Langley Research Center for his helpful suggestions and continuous 



Table of Contents 


1. Optimization of Structures to Satisfy a Flutter 

Velocity Constraint 1 

1.1 Velocity Gradient Search 2 

1.2 Projected Gradient Velocity Search ....... . 3 

1.3 Mass Gradient Search k 

1.4 Partial Derivatives of the Flutter Velocity ......... 5 

2. An Automated Procedure for Determining the Flutter Velocity ... 10 

Fi gure 1 11 

i ■ 

Appendix A - Published References and Abstracts 12 

Appendix S3 - "Optimization of Complex Structures to Satisfy 
Divergence Velocity Constraints" an unpublished 
Technical Note 15 





1. Optimization of Structures to Satisfy a Flutter Constraint 

In this section of the report a mathematical method for minimizing 
the structural mass of a lifting surface which is subject to a specified 
flutter velocity constraint will be presented. It will be assumed that 
structural parameters (cross-sectional areas, plate thicknesses, diameters 
squared, etc.) are selected in such a way that the total structural mass 
is a linear function of these parameters. The optimization procedure 
which will be presented is independent of the aerodynamic theory which 
might be selected. 
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1.1 Velocity Gradient Search 

This routine is employed when it is desired to increase the flutter 
velocity. The flutter velocity normal derivative is calculated at a point 
and parameters are varied so that a step is taken in the direction of 
maximum increase in velocity. The desired flutter velocity is reached 
after several successive steps along a gradient curve in an iterative 
fashion as from point A to point B in Fig. 1. 

The desired change in a design parameter during a velocity gradient 
search may be computed from the following relation which was derived in 

f 

Ref. 1. : 

n 2 

AP. = V . (V*-V)/E V . (1) 

' ’ j = l ,J 


where V* is the desired flutter velocity V is the current flutter 

velocity, V . = 3V/BP. and n is the total number of design parameters. 
,ii 

New values of the design parameters may be determined from the 
relation 

P? « P. + AP. (2) 

ill 


A new value of the reduced frequency k may be estimated from the 
relation 


n 

k* = k + Z k . AP. 

: _ i » ^ * 


(3) 


This new value of k may be used as a starting value for search for the 
true flutter velocity in a procedure for determining the flutter velocity. 
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1.2 Projected Gradient Velocity Search 

This routine is employed to find a relative maximum of the flutter 
velocity while the total mass of the structure is held constant. The 
parameters are varied in such a way that the search proceeds tangent to 
a constant mass hyperplane in the direction of the maximum rate of 
increase of the flutter velocity until a relative maximum is found which 
lies within the bounds of the parameter constraints as from B to C in 
Fig. 1. In Ref. 1 a procedure for searching for a maximum flutter 
velocity along a constant mass hyperplane was presented; however, the 
step size in the search was selected by trial, this resulted in a slow 
convergence to the maximum flutter velocity for that hyperplane. In 
Ref. 2 equations were presented for approximating the step size for 
moving from point B to C in Fig. 1. These relations are as follows: 


S 




/ 



n 

Z 

h=l 





where 



(V . + 6, m .)/20 
,J 1 ,J o 


n n 2 

e. - - I m . V ./ Z m 
1 j-i > J -j h=i ■' 


20 -■ 
o 



(V 2 .- + V .0. 
>J ,J 1 



1/2 



( 5 ) 

( 6 ) 

( 7 ) 


and where m is the total mass of the structure. New parameters may be 
computed from the relation 


P* 

J 


?. + (dP. /ds)S 


( 8 ) 
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1.3 Mass Gradient Search 

This routine is employed whenever it is desired to reduce the flutter 
velocity as from point C to point D in Fig. 1. A step is taken in the 
direction of the greatest rate of decrease in the structural mass and the 
process is repeated until the flutter velocity is less than or approximately 
equal to the desired value. 

The changes in the design parameters are computed from the following 
equation which was derived in Ref. 1: 


AP. 

i 


n 

m . (V*-V)/ Z V . m . 
» * ; — i > J » J 


(9) 


The above three search procedures were applied to the design of a box 
beam in Ref. 1 and 2. The results of this application indicate that 
practical application of the above method to a real aircraft structure is 
feasi ble. 
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1.4 Partial Derivatives of the Flutter Velocity 

Two of the search procedures which were previously described require 
the partial derivatives of the flutter velocity with respect to the design 
parameters. These partial derivatives may be found by considering the 
equations of oscillation of an aircraft structure which is in a state of 
neutral stability, i.e. 


[K - A. (M+A)]U. = 0 (10) 

where K, M and A are the stiffness matrix, mass matrix and aerodynamic 

f 

force matrix, respectively. U. is an eigenvector corresponding to the 
eigenvalue A.. The aerodynamic force matrix A is a function of the air 
dens i ty ? Mach number and reduced frequency k. 

Expressions for the first partial derivatives of the eigenvalues with 
respect to structural design parameters were derived in Ref. 1. These 
derivatives may be expressed in the form 


A. 


' >J 


= V' 


[K 


(M 


\i )] 


(ID 


where the letter in the subscript which follows the comma indicates a 
partial derivative with respect to a design parameter, i.e. 


A. . = 3A./3P. 

• »J i J 

The superscript T in equation (2) indicates the transpose of the left- 
hand eigenvector V where 


V. T [K - A. (M+A) ] = 0 

l I 


(12) 
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In the design process presented here it is convenient to constrain 
the derivatives of A. and k such that they are pure real. Then 
velocity gradient search and the projected gradient velocity search will 
move from one neutral stable state to another neutral stable state for 
the mode which determines the flutter velocity. The derivatives for 
this case are given in Ref. 1 and Ref. 2 and are repeated here. 



X. . — R. — R_ A./A„ 

i , j 1 2 12 

(13) 


k . = - A./A 

> J ' *■ 

(14) 

where 


! 

R 1 + 

A, I = vT [K . - A.M .]U. 
1 m J ,J i ,J 1 

(15) 

*2 + 

. , , ,,T 3A .. 

A 0 1* = - A. V. -rr U. 

2 m i i 3k i 

(16) 

1 - 

<-o ,/2 


m 


The f i rst partial 

derivative of the flutter velocity 

i s g i ven 

by the relation 



V fJ 

= - boj. k ./k^ + bA. ./2kw. 
• ,J i »J i 

(17) 


where |/ 2 

CO. = A. 

\ I 

The second derivatives of the eigenvalues of equation (10) and the 
reduced frequency when the derivatives are constrained to be pure real are 
derived in Ref. 2 and fray be expressed in the forms 

x i ,jh ° R 3 ' *2 V A 2 (l8) 
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and 


09 ) 


where 


\jh ‘ " A 3 /A 2 


R, + A, I = V! < 
3 3 m i 


K .. - X. . (M . + A , ) 
,jh ' *J ,h ,h' 


- X. . (M . + A .) - X. [M ..+ (3 2 A/s|)/k .k . ]>U. (20) 

I ,J ,J < ,jh ,J , h ( i 


+ vT F. . U. . + vT F. . U. . 
i i,J i,h i i,h i ,j 


A . = 
»J 


A 

>jh 


3A . 

9k \j 


5 2 A 

3k 2 


k . k , + 


3A 

8k 


k 

,jh 


( 21 ) 
. ( 22 ) 


v: f. .o. . 
i i ,j i ,h 


T ^ T 

= V [K .“X (M .+A .)]E V [K ■ -X (M ,+A ) ] U . U / (X -X, ) 

I ,J l ,j ,J * ,n ? ,h ,h i 1 f 1 


(23) 


1^1 


The second derivatives of the flutter velocity with respect to the 
design parameters may be found by differentiating equation (17) with 
respect to a design parameter, then 


V = (b/2kw.) [X. .. - X X. , / (2X. ) 
ijjh i i,j i > n i 


-(X. .k + X. k .)/k 
i ,j ,h i ,h ,j 


+ (2k .k . / k - k .. )2 X./k] 
,j ,h »jh i 


(2A) 


Alternate methods for finding the unconstrained derivatives of eigen- 
values and eigenvectors were developed in Ref* ^ and 5* The "algebraic 
method" of Ref* 5 rnay be used to find the first derivatives of the eigen~ 
values and eigenvector of the non-self-adjoint system of e i genequat I on 


[K - A * (M+A)]U. = 0 


( 10 ) 
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with respect to an independent variable P. 
equations 



by solving the set linear 


( 25 ) 


where 



K - A. (M+A) j 


0 1 


- (m+a) u . 


( 26 ) 


and 

' D = 

and where it is required that ulu. = 1. 

i i 

The second partial derivatives of the eigenvalues and eigenvectors 
may be found by differentiating equation (25) with respect to the inde- 
pendent parameter P^ , then 


A. (M . + A . ) - K . 


( 27 ) 



( 28 ) 


If A. . , U. . and U. , are evaluated from equation (25) then 
* * J * > J 1 > 

U. ♦, and A. may be found by solving equation (28). 

1 , j n 1 > J n 

The above process may be continued for any number of higher order 
der i vat i ves . 

It should be noted that the derivatives found by using the algebraic 
method are not constrained, and they may not be substituted into 
equations (i7) and (2A) to find the first and second derivatives of the 
flutter velocity; however if the reduced frequency is held constant in 
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equations (25) through ( 28 ) and in the relations 


V f = b to./ k = b \\ /2 /k 


then 


\l f . = b A. ./2ku. 
f .J i,j i 


and 


v fjh = b A i >jh /2k "i - b A i,j A i >h ^K 


(29) 

(30) 

(31) 


The derivatives of V f from equation (30) and (31) will (in the 
general case) be complex. 
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2. An Automated Procedure for Determining the Flutter Velocity 

Efficient optimal design programs for aircraft structures which are 
subject to constraints on the flutter velocity require a rapid and 
automatic method for evaluating the flutter velocity. In Ref. 3 a 
computationally efficient method for finding the flutter velocity is 
presented. The method utilizes derivatives of the eigenvalues with 
respect to the reduced frequency in a curve fitting scheme for finding 
the critical roots of the flutter equation. The method Is unaffected 
by the coalescence of any of the eigenvalues. 

In Ref. 3 the derivatives of the eigenvalues were found by use of 
the methods of Ref. 1 and 2; however, the "algebraic" method of Ref. 5 
could be more efficiently employed to find these derivatives. 
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Abstract 

Equations for finding the partial derivatives of the flutter 
velocity of an aircraft structure with respect to structural 
parameters are derived. A numerical procedure is developed for 
determining the values of the structural parameters such that a 
specified flutter velocity constraint is satisfied and the 
structural mass is a relative minimum. A search procedure is 
presented which utilizes two gradient search methods and a 
gradient projection method. The procedure is applied to the 
design of a box beam. 
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Mo 1 . 10, No. 12, December 1972, pp. 1569-1572. 

Abstract 

Equations for the second partial derivatives of the eigenvalues of 
the flutter equation along with the equations for finding the second 
partial derivatives of the flutter velocity of an aircraft structure 
with respect to the structural parameters are derived. These partial 
derivatives are used to develop expressions for the step size in a 
projected gradient search along a constant mass hyperlane. A 
projected gradient search along with a gradient mass and a gradient 
velocity search is used to minimize the mass of a box beam which 
supports a lifting surface. 

3. Rudisili, C.S. and Cooper, J.L., "An Automated Procedure for Determining 
the Flutter Veloc ity ‘ Journal of Aircraft, Vol. 10, No. 7, July 1973, 
pp. 442-444. 


• . Abstract , .. 

A computationally efficient method for finding the flutter velocity 
of a structure is presented. The method utilizes derivatives of the 
eigenvalues with respect to the reduced frequency in a curve fitting 
scheme to find the critical roots of the flutter characteristic equation. 
The method is unaffected by the coalescence of any of the eigenvalues. 



4. Rudisill, C.S., “Derivatives of Eigenvalues and Eigenvectors for a 
Genera] Matrix," AIAA Journal , Vol. 12, No. 5> May 1974, pp- 721-722. 

Abstract 

A method is developed for finding all orders of the derivatives 
of the eigenvalues and eigenvectors of a non-self-adjoint system of 
algebraic eigenvalue equations. The method does not require a 
complete ei genana 1 ys i s of the algebraic eigenvalue problem. The 
method requires only the eigenvalue and eigenvector whose derivatives 
are sought and the method also requires the corresponding left-hand 
eigenvector. 

5* Rudisill, C.S, and Chu, Yee-yeen, “Numerical Methods for Evaluating 
the Derivatives of Eigenvalues and Eigenvectors," Accepted for 
publication in the.AIAA Journal. Publication date will probably be 
in the late spring of 1975. , 


Abst ract 

Two numerical methods for computing the derivatives of eigenvalues 
and eigenvectors are developed. The first method is an iteration 
method for finding the first partial derivative of the eigenvalues 
and eigenvectors of a self-adjoint system of algebraic eigenvalue 
equations. The iteration method will also find the first partial 
derivative of the largest eigenvalue and its corresponding eigen- 
vector of a non-self-adjoint system, but the method cannot be used 
to find the derivatives of the remaining eigenvectors. 

The second method will find all orders of the derivatives of the 
eigenvalues and eigenvectors of a non-self-adjoint system of algebraic 
eigenvalue equations. The method does not require a complete eigen- 
analysis of the algebraic eigenvalue problem. The method requires 
only the eigenvalue and eigenvector whose derivatives are sought. 
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TO SATISFY DIVERGENCE VELOCITY CONSTRAINT' 
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ABSTRACT 

I 

Analytical expression for evaluating the partial derivatives of 
the torsional divergence velocity of an aircraft structure with respect 
to design variables is derived. An optimization procedure to satisfy 
a specified divergence velocity is i 1 1 ustrated, usi ng gradient methods 
and finite-element representation, for a box beam with the lower and 
upper values specified for the design variables. It is shown that 
there is a possibility of serious designer judgement error due to 
inefficient performance of optimization methods against multiple 
constraints. A "dimension reduction technique" is proposed to help 
in such situations. 
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1. INTRODUCTION 


In the recent past, there have been several publications dealing 
with optimization with respect to a static aeroelastic constraint on 
the torsional divergence velocity [1,2, 3, 4] • McIntosh and Eastep [2] 
presented a calculus of variations formulation for a tapered cantilevered 
wing with torsional stiffness dominated by contributions from a thin 
outer skin. McIntosh and Veisshhaar [3] > and Armand and Vitte [4] used 
transition-matrix approach of optimal control theory. McIntosh and 
Weisshhaar [3] concluded that early optimism concerning the use of 
transition-matrix approach must be tempered somewhat, and it may in 
the long run prove most useful to adopt more sophisticated steepest 
ascent or gradient methods. Further, it appears that even for the 
gradient methods to be computationally attractive for practical problems, 
v/1 thout excessive penalty in the computer time used to arrive at an 
optimum solution, closed form analytical expressions should be derived 
for the partial derivat!ve.s involved so that their computation does not 
involve individual perturbation of each design variable and corresponding 
evaluation of the objective function. This would be a very significant 
consideration for problems with a large number of design variables, 
and/or where evaluation of the objective function is expensive in terms 
of computer time. In the present paper, the authors' aim Is, therefore, 


to develop a. closed form expression for. the partial derivatives pf. ... . , 

torsional divergence velocity with respect to design variables, and to 
present a practical method for optimum weight design of an aircraft wing 
structure subject to torsional divergence velocity constraint. 


Numbers in square brackets indicate references listed at the 
end of paper. 
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H. PARTIAL DERIVATIVES OF DIVERGENCE VELOCITY 

The governing characteristic equation for torsional divergence 
of an aircraft wing structure in a static neutral state can be 
expressed as, (see Ref. (5 )) 

(tKj - A [A]) {U} =0 (1) 

where 

[K] = torsional stiffness matrix, function of design 

variables P., symmetric matrix. 

[A] = torsional air-force matrix, function of air 

density and wing geometry, constant and real . 

s 

A = eigenvalue, equal to the divergence velocity 
squared. 

{U} = angular displacement of the wing 

To define an associated row vector {V}' of the eigenvector (U), 
consider 

([K] - A [A] ' ) {V} = 0 (2a) 

Taking the transpose of equation (2a) and using the symmetry 
property of [K] yields 

{Vj'UK) - A[A]) = 0 (2) 

Differentiation of the character i st ic equation (l) with respect 
to a design variable P. yields 

) : ;u) * (IK] ..-.X (ah o. .. 

i i i 

Premultiplying equation (3) by {V} and simplifying by substituting 
equation (2), results in the equation 
(V) ' (U) 

3A = ■ (**) 

3P i (V)'(A] (U) 


3 



A similar expression can be derived using the flexibility 
matrix instead of the stiffness matrix, in which case the characteristic 
equation is given by 

( [C] [A] - 1A[I]) {U>- 0 (5) 

where [C] is the inverse of [K]. 

From equation (5) it can be shown that 

x 2 <u) (6) 

3X - JJ ' 

3P « {V}'{U} 

The partial derivatives of the divergence velocity can be 
quickly computed from either equation {*0 or (6). the use of these 
derivatives v/ill be discussed in section III. 



HI. OPTIMIZATION PROCEDURE 


• Two genera} approaches to optimization have e^volved thus far: 
maximizing a certain eigenvalue for fixed mass, or minimizing the mass 
for a fixed eigenvalue [2]. The authors use the first approach, which 
is equivalent to projected gradient method wi th eigenvalue as the 
objective function, in conjunction with gradient mass and gradient 
veloc i ty methods . This procedure ?s described in detail in the 
authors' earlier paper [6]. The optimum design with divergence 
velocity constraint will be illustrated for a three bay box beam 

t 

.(See Figure 1) with twelve design variables having upper and lower 
constraints. 

The optimization program is independent of the formulation used 
for the air-force matrix. However, significant simplification is 
obtained when the air-force matrix [A] is symmetric, in which 
case equation (4) is simplified to 

iv}'[ f£- r { u> 

B - \ (7) 

3P l {U}‘[AJ {U} 

This simplification would not be achieved using the flexibility 
matrix. Thus there is substantial computational advantage for the case 
of symmetric air-force matrix, if the stiffness approach is used. 

For the numerical example presented in this paper, a diagonal air-force 
matrix was used [5]. 

Table 1 lists the results obtained from the optimization program 
for the box-beam of Figure 1. Two sets of values for the design 


s 



variables were used as Initial input to the optimization program. 

The mass-optimization was subject to a torsional divergence velocity 
constraint of 600 ft. /sec. with upper and lower constraints specified 
for all the 12 design variables. These constraints were same for 
both the cases of initial design values. 

For the case 1, an optimum design was very quickly arrived at 
m two design cycles. A design cycle here defines a step taken in the 
mult i -dimensional space of design variables during an execution of 
projected gradient, mass gradient or velocity gradient search; each 
design cycle involves computation of stiffness matrix and evaluation 
of the eigenvalue problem in addition to computation of derivatives. 

The final values for P 2 , P 3 and P ]Q seem to be numerically same as 
the corresponding initial values. Actually the final values are 
slightly higher than the initial values, but due to a very small 
numerical difference this is not apparent from the number of digits 
listed in the table. Therefore, the optimum arrived at for case 1 , 

dorvs t y^.' * °£ k e a ^ ree °Pt imum ; parameters <xre, net a.ga<r>st any 
For the case 2, the lower constraints for P^, p g and P 
slewed down the optimization procedure considerably, and ittook 4? 
design cycles to reach a mass approximately Ii0% higher and a divergence 
velocity approximately O.W higher than for the previous optimum attained. 

This suggests the possibility of serious designer judgement error 
duo to inefficient performance of optimization methods against 
multiple constraints. To circumvent this, the authors suggest the 
use of "dimension reduction technique". The proposed technique would 
monitor the partial derivatives of the relevant eigenvalue with respect 
.to the design parameters which are against the constraints, and from 



these determine the design parameter which would tend to violate the 
constraints if a step in the desired direction were taken. Such 
parameters would be then held constant for the next step, thus 
effectively reducing the dimension of the design parameter space 
for this step. At the new point thus reached, a new design cycle 
would begin and partial derivatives of the relevant eigenvalue 
with respect to all the design variables would be computed, and 
the above process repeated. It is expected that such a technique 
would reduce the number of design cycles required for problems 
where constraints are encountered. 



CONCLUSIONS 


The closed form analytical expression derived for the partial 
derivatives of divergence velocity with respect to a design 
parameter is useful in the gradient type search procedures. The 
example solution illustrates its use. It seems that where lower 
constraints are specified for the design variables, a logical 
initial point should preferably include applicable lower constraints. 
This view is substantiated by authors experience that gradient 
velocity procedure operating In the neighborhood of lower . 

constraints is very fast and effective in increasing the velocity 
to the desired value, since in increasing the velocity It^ould 
usually tend to move away from the constraints. For the case where the 
constraints slow down the gradient methods, the dimension reduction 
technique may improve their performance. The authors hope to explore 
the potential of this method in their future work. 
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TABLE 1 Initial and final Parameters with Constraints 
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— 

of Longitudinals 
(So. in.) 

Front £ Back Web Thickness 
(x 10 1 in . ) 

. 

Top £ Bottom Web Thickness 
(x 10 ^ in.) 

Rib Thickness 
(x 10 _1 in.) 

Tota i 
.Mass , 
,( S 1 u g 5 , 

1 

!D i verger.ce 
jVe loc ; ty 

! (ft ./ sec . ) 




Bay 2 
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P 3 

Bay 1 
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Bay 1 

P 7 

Bay 2 
P 8 
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0.0666 

0.0666 
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0.S2 

555.23 

— a 

C f? 

— > 
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lA'Cst 

O 
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2.000 

2.000 
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‘ 
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■ 

■ 

0.400 : 

■ 

0.400 

0.400 

0.400 



1370.41* 

w 

P 

■ 

0.33264 

0.33264 

0.33264 

0.13332 

0.13332; 

0.13332 

0.0666 ' 

• 

0.0666 

0 * 0666 

0.0666 

0.0666 

0 . 0666 

0.52 . 

555.23 

0 
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k- 

o 

CL 
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Z3 

8.0C64 
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8.0064 
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3.204 . : 
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■ 
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— 

22.15 

.2733.71* 

tn 

— O 
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iA 

Q ^ 
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0.33336 


0.33264 

0.1362 
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0,133368 

0.081228 

0.082284 

0.072672 

0 . 0666 

0.066756 

0.067008 

0.95 

595. SS 

1 

fU 13 

L- > 

O 

LO 

ro cm 

0.3865 

0.3865 

< j 

O.3865 

0.35184 

0.35184 

0.3184 

0.0666 

0 . 0666 

0.0666 

0.3066 

0.3066 

0.3066 

1.30 

; 602.65 


* Beyond the range of aerodynamic theory used, 
a) Number of design cycles: Case 1-2 

Case 2-47 
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